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^—^ i Abstract. The propagation of a gravitational wave (GW) through a magnetized plasma 
is considered. In particular, we study the excitation of fast magnetosonic waves (MSW) by 



o 
S3 

Oh 
i 

o 



a gravitational wave, using the linearized general-relativistic hydromagnetic equations. 
We derive the dispersion relation for the plasma, treating the gravitational wave as 
a perturbation in a Minkowski background space-time. We show that the presence of 
gravitational waves will drive magnetosonic waves in the plasma and discuss the potential 
^ astrophysical implications. 
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1. Introduction 

It is well known that gravitational waves (GWs) can interact with electromagnetic waves. 
The scattering of electromagnetic waves by time dependent GWs has also been considered 
by several authors (DeWitt & Breheme 1960; Cooperstock 1968; Chestrs, 1973; Denisov 
1978; Grischuk & Polnarev 1980; Macdonald & Thome 1982; Machedo & Nelson 1983; 
Demianski 1983; Marklund et al. 2000). It has been shown that the coupling is mainly 
due to the generation of electric currents inside the plasma by the perturbation of the 
charged particle trajectories by the GW. Several studies used the kinetic theory formula- 
tion to study the coupling of GWs with electromagnetic waves in a magnetized vacuum 
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(Gertsenshtein 1962; Zeldovich 1974) or near a charged black hole (Gerlach 1974). The 
first attempt to study the propagation of GWs in a magnetized plasma using the kinetic 
theory was by Macedo and Nelson (1983). They studied the coupling of GWs propagating 
parallel and perpendicular to the external magnetic field with high-frequency electromag- 
netic modes. The coupling of the gravitational with the high-frequency electromagnetic 
modes propagating along the magnetic field was weak. This is due to the fact that the 
GW propagating along the direction of the undisturbed magnetic field does not generate 
currents in the plasma. The coupling is substantially higher for ordinary and extraor- 
dinary waves propagating perpendicular to the magnetic field but the energy transfer 
of gravitational energy to electromagnetic energy in the interstellar space is extremely 
weak. 

Brodin & Marklund (1999) discuss the parametric excitation of plasma waves in the 
presence of gravitational radiation. They assumed that the GW was a monochromatic 
wave and estimated, using Einstein-Maxwell- Vlasov's equations, the growth of Langmuir 
and electromagnetic waves in a cold unmagnetized plasma. Marklund et al. (2000) used 
the 1+3 orthonormal frame formalism (Ellis & van Elst, 1998) to study the linear inter- 
action of a GW inside a inhomogeneous magnetized plasma. They showed that, close to a 
strong gravitational source, the electric fields associated with the electrostatic waves can 
reach appreciable levels. The direct escape of plasma radiation is not possible without 
a non- linear conversion (see Brodin et al. 2000). Recently Servin et al. (2000) discuss 
the non-linear parametric excitation of Alfven waves by a strong monochromatic GW. A 
nonlinear study of shock waves in a highly magnetized plasma, driven by GWs has been 
presented by Ignat'ev (1997). 

The non-linear interaction of a monochromatic GW with charged particles and their 
cyclotron damping inside a plasma has also been analyzed in detail (see Esposito 1971; 
Polnarev 1972; Asseo et al. 1976; Papadopoulos & Esposito 1985; Varvoglis and Pa- 
padopoulos 1992; Servin et al. 2001) 

Papadopoulos and Esposito (1982) derived the exact equations for a finite amplitude 
MHD wave propagating inside the plasma, within the framework of general relativity. 

In this article we study, for the first time, the excitation of magnetosonic waves from a 
GW propagating inside a uniform, strongly magnetized plasma (where the Alfven speed 
approaches the speed of light). In section 2 wc outline the basic equations, in section 3 we 
derive the general dispersion relation for the coupling of the GWs with the MHD waves 
and in section 4 the dispersion relation for the coupling of the gravitational with ambient 
magnetosonic waves. In section 5 we discuss our results and their potential astrophysical 
implications, while section 6 summarizes our conclusions. 
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2. Basic equations 

The equations for finite amplitude wave propagation in an ideal hydromagnetic plasma 
in the general relativistic framework are (see Papadopoulos & Esposito, 1982) 



(e - ^-);a b u a u b = h ab { P +^-), ab + 2{H 2 Q)- 

9P> 2 

- (H a H b ), ab + 2x(-f- + a 2 -u 2 

- u a u a ) + -(p + 3p 

+ H 2 ) + 2u a {H a H b ). b + (H 2 ). a u a . (1) 

Maxwell's equations take on the compact form 

H a = (a a b + cot - h a b Q)H b + -^- P[b H b u a , (2) 

S e + p 

and the equation of motion is 

H 2 
xu a + xu a + xQu a + (p+ ^-);h.9 afc - (H a H b ). b = (3) 

where u a is the fluid velocity, e is the energy density, p the rest mass density, p the 
pressure of the fluid, H a = (0, H 1 , H 2 , H 3 ) is the prevailing magnetic field, ii a = u a c u c , 
x = e + p + H 2 , h ab = g ab + u a u b , Oij is the shear, w tJ is an antisymmetric tensor, 
corresponding to the rotation of the fluid and O is the expansion of the fluid. In the 
equations we have assumed c — 1, k = 8irG/c 4 = 1 and we have put H = B —> \/4ttH 
(from Gaussian units), where k is the coupling constant in Einstein's field equations. 

We perturb the equations (l)-(2) and keep only linear terms in the perturbed quan- 
tities (see Eqs. (47)- (49) in Papadopoulos and Esposito, 1982), and use the metric given 
by g^ v = V^u + hp V or 

ds 2 = -c 2 dt 2 + (1 + h)dx 2 + (1 - h)dy 2 + dz 2 (4) 

where h = h(t, z) represents the perturbation of the flat space by the GW. 

In the following section we will search for the coupling of the normal modes of the 
magnetized plasma with the GW. It is important to notice that both gravitational and 
MHD waves follow similar dispersion relations, in the asymptotic limit when the rela- 
tivistic Alfven speed approaches the speed of light (see also Landau and Lifshitz, 1975). 

3. The Perturbed Equations 

In this section we will perturb Eqs.(l)-(2) assuming that Sg ab ^ 0. Perturbing Eq.(l) we 
find 
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+ (H 2 ), a (5u a ), c u c 
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The perturbed Maxwell's equations are 

5H b + ST ol H l = -Sa b a H a + ^H b SQ - -^—Sp tC H c (6) 

and from Eq.(3) we find the perturbed equations of motion 



- (e + p + H 2 )Su% = h a [H b 5H% + 5T a bl H b H l 

+ H a 5H b b + ST b bm H m H a ] - 5h c % - h cb S P : b (7) 

where p* = p+ ^-. The non-zero Christoffcl symbols T" , for the metric (4), are 

estimated easilv V° - -T° - _^£ r° — T° - /l -° T 1 - T 1 - h ' 3 
estimated easily i n — i 22 — 2 , i 01 — i 10 — 2(l+fc) ' 13 — L 31 — 2(l+h) ' 

T 2 2 = r| = -g^, T 2 3 = r| 2 = -5^, T 3 U = T\ 2 = *£. We assume that 
u a = (-1,0,0,0), p = / o =const., H a = (0,H 1 ,0,H 3 )=const., 8H a = (0, SH 1 , SH 2 , 0), 
(5w a = (<5w°, <5w 1 ,(5'u 2 ,(5'U 3 ) and k a = (0,0,0, A:) and a plane wave behaviour for the per- 
turbed quantities Se, Sp, Su a and SH a ~ ^(kz-nt) an( j ^e GW h ~ e *(M-™ 3 *). These 
simplifications yield a general dispersion relation for the coupling of GW with MHD 
waves. 



Se{D (-n 2 - 2n 2 u z A ) + (kc s ) 2 D Q 
±(l + c 2 s )(e + 3p + H 2 )D 



D. Papadopoulos et al....: Magnetosonic waves driven by Gravitational waves 



- Ul + 3c 2 s )(e+p + H 2 )D Q 



1„2\ 



[n 2 + k 2 - 2n 2 u\ - 2(e + 2p + iJ 2 )](AiiJ 1 + X 2 H 

[n 2 - k 2 - 2n 2 u 2 A - 2(e + 2p + H 2 )]( Ml ff 3 + /i 2 ff 3 c 2 )} 

hiDoiH 1 ) 2 ^ + 2 P + H 2 - k 2 ) - XzH 1 - ^H 3 } 



(8) 
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Above, 



fi/(i 



is a relativistic generalization of the classical Alfven speed 



va = H 2 /{Aite/c 2 ) and c s = y/Sp/Se is the speed of sound. In passing we note that the 
non-relativistic Alfven speed va generalizes to the relativistic Alfven speed ua if one does 
not neglect the displacement current in Maxwell's equation, in MHD theory. 
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4. Magnetosonic Waves 

We are interested in investigating the dispersion relation for waves that correspond to the 
magnetosonic mode (H 1 ^ 0,H 3 = 0). This is motivated by the fact that the coupling 
of MHD waves with GWs is proportional to (H 1 ) 2 . Thus, gravitational waves in the 
z-direction cannot couple to pure Alfven waves (H 1 = 0, H 3 ^ 0). The general relations 
between the perturbed quantities, derived above, simplify when we assume that H 3 = 
(setting H 1 = H from now on) 



H 2 1 

-n 2 5e - 2n 2 Se + k 2 Sp - -(Se + 5p)(e + 3p + H 2 ) 

x 2 

- ^(Se + 3Sp)(e+p + H 2 ) + (HSH ) [n 2 + k 2 

H 2 

- 2n 2 2(e + 2p + H 2 )} 

x 

= h(H 2 )[(e + 2p + H 2 )-k 2 ] (9) 

where, from Maxwell's equations (6) we obtain 

{H5H) = -^[--fpL\. (10) 

1 — u A x n 

From here on, we will write our equations in ordinary units again. We combine Eqs.(9)- 
(10) to obtain 



5e{-n 2 + (kc s ) 2 + k 2 u\ - (kc a ) 



2 U A 



H 2 

c 2 ' c- c- 

H 2 
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- 4G^(47rp+— )} 
c z c 

= h^{^[n 2 (l-24) + c 2 k 2 -AG(^p+^)} 



+ 2(l-^)(4nGp+^-k 2 g c 2 )} (11) 

where H 2 = H a H a = H\H l and we have set p = and p = e/c 2 . 

We retain the dominant terms that yield a pure magnetosonic wave in a strongly 
magnetized, cold and tenuous plasma (c s << Ua) in the Newtonian limit. This is cquiv- 

4/ 
obtains 



alent to setting 1/c 2 — > (except for the u 2 A /c 2 terms) and Gp — > 0, c s — > and one 



(-n 2 + k 2 u 2 A )Se = h^-i^- 
4n n 



>r(l-2f 
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k 2 c 2 -2(1-^ U V . (12) 
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In the absence of a GW (h = 0), the dispersion relation ( fl2| ) becomes 

n 2 r ~ fc 2 ^- (13) 

If one assumes that the frequency of the driving GW coincides with the frequency of the 
produced MSW, i.e. n = n g = k g c, then the dispersion relation becomes 

(-n 2 g + k 2 u 2 A )5e = h^-(-n 2 g + k 2 c 2 ), (14) 

or 

I? " *l?) * = "£(-"» + fcV >- <15 > 

i.e. the GW acts as an external parametric oscillator, driving the MSW. Notice also that, 
in the absence of a plasma one obtains the dispersion relation for GWs propagating in 
the vacuum. The two wave numbers differ by an amount Afc, i.e. k = k g + Afc, which, in 
general, is complex, Afc = Afc/j + iAfc/. If we denote 
f h\ H 2 

(where a subscript denotes amplitude), then (|14|) yields a system of two equations, for 
Ak R and Afc/. Dividing the two resulting equations, one obtains 
-n 2 g + [k 2 + 2k g Ak R + Ak%] u\ A cos + B sin c\> 



2k g Ak]u\ A sin + B cos (j> ' 

where 



(17) 



A = [2k g Ak R + Ak 2 R ] c 2 , (18) 

B = 2k g Ak lC 2 . (19) 

and 

<p = Ak R z. (20) 

When | Afc/ 1 << |Afc#| =>■ |S/j4| << 1 and since ua — c in a strongly magnetized, 
tenuous plasma, we arrive at an expression for Afc/ 

Afc/~-^L(l-^)tanAfc/^. (21) 

The real part of the wavenumber difference is approximately 

Ak R = k- k g = k g (^--lj, (22) 



thus, the two waves can linearly interact coherently over a maximum distance of 
A„ 



27r(l — Mq/c) 
in terms of the GW wavelength X g . 



(23) 
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For the plasma around a magnetized compact object, B 2 /e ~ 10 2 — 10 4 (where e 
is dominated by the kinetic energy of the plasma particles). For such a plasma, the 
interaction length can be up to several thousand times larger than X g , A more stringent 
limit on the interaction length will probably be set but factors such as the non-uniformity 
of the magnetic field and the plasma density with distance from the source. 

Assuming also a small imaginary part for the frequency and substituting n = n r + ij 
in Eq.(14), developing — n 2 + k 2 u 2 A = (n — n r + *7)(— 2n r ), and using the Plemelj formula 
( _In+ i x + *7) _1 — P{^-/ x ) ~ i^5{x)) we finally arrive at the spatial dependence of the 
amplitude of the energy density 

Se » U{n - n g )h H 2 c^ke- l(Ak)z , (24) 

where Se = 5e e^ nt ' kz \ SH = SH e l ( nt ~ kz *> and h = h e l ^ n ^- k ^\ 

The magnetic field amplitude is estimated with the use of the perturbed Maxwell's 
equations (6) and with the aid of Eq.(lO) 

SH m l -8{n - n g )h HcAke- lAkz . (25) 

We note that, the amplitude of the fast magnetosonic wave is proportional to the GW 
amplitude. In a tenuous and hot plasma, ua — ► c, so that the mismatch A/c « k and our 
treatment applies. As a result the amplitude of the MSW grows after having traveled a 
distance equal to a half beat between the wave numbers, and would decline again further 
out in the absence of dissipation. In a relativistic electron positron plasma the latter is, 
however, unlikely and some observational effects can be expected. 

5. Discussion 

5.1. Interpretation 

We have found that, to first order, a GW, propagating along the magnetic field in a 
uniform and ideal plasma, does not couple with an Alfven wave (AW) whereas a GW, 
propagating across the magnetic field, does couple with a MSW. This can be easily be 
understood as follows. 

In Fig. |] we have sketched the effect of a large-amplitude GW traveling along an, 
initially uniform, magnetic field embedded in an ideal plasma. At first sight it is sur- 
prising that the GW would not be able to excite an AW as the perturbed magnetic 
field shows clear distortions which in an ordinary magnetoplasma would be expected to 
generate AWs. However, the point is that the energy of an AW with amplitude 5H l 
is proportional to (SH 1 ) 2 whereas, for a MSW, the magnetic part of the wave energy 
density is proportional to H^^SH 1 . As a result, the energy in the AW is of higher order 
in the perturbation as compared to the magnetosonic case. Therefore, our first order 
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calculation should not lead to the excitation of AWs, as indeed we have found. On the 
other hand, we do expect such a coupling to show up in a higher-order calculation, and 
such a coupling will be relevant in case the GW is non- linear. 

In Fig. H we have sketched the effect of a GW traveling at right angles to an initially 
uniform magnetic field frozen into an ideal plasma. Clearly, the field is periodically com- 
pressed by the GW although the differential motions induced by the GW are divergence- 
free! At the same time, the ambient uniform magnetic field remains uniform and straight 
and does not become bent. As a result, no AW excitation is expected to occur by the 
GW even to higher order. 

Further, our result that the coupling between a GW and MSWs vanishes as soon as the 
dispersion relations for both wave modes coincide exactly, can be understood as follows. 
An effective transfer of energy from one wave mode to the other either requires that there 
is some dissipation (Weber 1960; Esposito 1971) or, alternatively, that the signal speeds 
of both modes differ. If the plasma is both ideal and the group speeds of the wave modes 
coincide, there is no time lag during which some of the perturbed kinetic/magnetic plasma 
energy established by the GW can be fed into a magnetosonic wave. In a very general 
sense, the changes brought about by the GW in the magnetoplasma are completely 
adiabatic and the periodic changes in the plasma completely reversible. All free energy 
in the magnetic field reverts back into the GW as it leaves the plasma again. 

5. 2. Applications and perspective 

A number of papers exist which discuss the interaction of MSWs with a plasma, e.g. 
coherent MSW-particle interactions, cascading of long wave-length MSWs to short wave- 
lengths and subsequent particle acceleration from turbulent diffusion (sec Dcrmer at 
al. 1996), and finally the interaction of a relativistic magnetosonic shock wave with a 
positron-electron plasma (Hoshino et al. 1992). Thus, if a magnetosonic wave is driven 
by a GW, then it has the potential of producing additional interesting phenomena. A 
more detailed understanding of the linear and nonlinear evolution of the excited MSWs 
under the conditions of the particular GW source is needed. Clearly, it is beyond the 
scope of the present study to model the acceleration and radiation properties of a GW- 
driven MSW in the vicinity of a GW-source but we want to mention briefly two specific 
applications which have been the motivation for us to undertake this study: GWs excited 
by star quakes or rotational instabilities in neutron stars (see Andersson, Kokkotas & 
Schutz 1999) or magnetars, and GWs excited in merging neutron star binaries and their 
possible effects on gamma ray bursts (GRBs). 

In particular, we expect that the driving of MSW by gravitational waves can become 
interesting near strongly-magnetized compact sources of high-frequency GWs. The non- 
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Fig. 1. Sketch of the effect of a large-amplitude GW traveling parallel to a uniform 
magnetic field frozen into an ambient ideal plasma, k \\ Bo at an arbitrary time. Here the 
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Fig. 2. Sketch of the effect of a GW traveling perpendicular to a uniform magnetic field 
frozen into an ideal plasma. The tidal motion of the GW is indicated with curved dashed 
lines. The drawn and, respectively, dashed arrowed lines are the magnetic field lines half 
a cycle apart. Their individual relative displacements are shown by the bold sections 
along the tidal curves. Note that the field lines remain straight under the action of the 
GW. 



radial modes of relativistic stars produce GWs with frequencies of the order of kHz. 
These can be excited in several ways: in the core-collapse leading to the formation of a 
neutron star, several modes of pulsation can be excited to large amplitudes. In newly- 
born neutron stars the quadrupole /- and r-modes can become unstable to the emission 
of gravitational radiation, due to the CFS instability (see Stcrgioulas 1998, Andersson 
& Kokkotas 2001 for recent reviews). In the merger of two neutron stars, if a black hole 
is not formed promptly, the resulting merged object will emit large amplitude GWs, 
due to the excitation of it's nonradial modes (see Shibata and Uryu, 2001). If strongly 
magnetized neutron stars (magnetars) exists, then the coupling of their secular spin- 
down with their very large magnetic field can lead to crust-fracturing and excitation of 
the quadrupole torsional mode in the crust of the star (see Duncan 1998). In all the 
above cases, high-frequency gravitational waves travel through the magnetosphcre of the 
compact object. 

In addition to GWs emitted by normal-mode excitation in neutron stars, the main 
frequency of GWs produced in the late stages of a neutron star binary inspiral is several 
hundred Hz to a few kHz (see e.g. Duez, Baumgarte and Shapiro 2000). In such a case, 
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however, the geometry of the magnetic field changes rapidly and this has to be taken 
into account. 

6. Summary and Conclusion 

In this article we study for the first time the coupling of magnetosonic waves with GWs. 
Using the linearized general-relativistic hydromagnetic equations, we found the general 
dispersion relation for coupling of GW with the MHD waves. We showed that in the 
linear approximation only the fast magnetosonic waves will be coupled with the GW. 
We analyzed the linear dispersion relation for the GW-MSW coupling, and found the 
relation between the amplitudes of the two waves, when they are in resonance. 

In carrying out this study we have made several simplification: (1) we treat only one 
polarization for the GW (see Eq. (4)), (2) the GW was treated as an external driver for 
the plasma, (3) we ignore non-linear terms and obliquely propagating waves. The above 
points, as well as the growth rate of driven MSW waves for specific GW sources, will be 
discussed in future work. 
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